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ABSTRACT 


We have applied the first Born approximation (FBA) to 

the calculation of single electron loss cross sections for various 

ions and atoms containing from one to seven electrons. Screened 

hydrogenic wave functions were used for the states of the electron 

ejected from the projectile and Hartree-Fock elastic and incoherent 

scattering factors were used to describe the target. The effect of 

the target atom on the scaling of projectile ionization cross 

sections- with respect to the projectile nuclear charge was explored 

in the case of hydrogen-like ions. We also have examined scaling 

of the cross section with respect to the target nuclear charge for 

+25 

electron loss by Fe in collision with neutral atoms ranging from 
H to Fe. These results were compared to those of the binary 
encounter approximation (BEA) and to the FBA for the case of 
ionization by completely stripped target ions. We have also 
calculated electron loss cross sections for the ions 0 (i=3-7) 

and N +1 (i=0-6) in collision with He targets in the energy range 
of ^0.1 to 100 MeV/nucleon. We found these results to be in excellent 
agreement with the available data near the peak of the ionization 
cross section. 
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I. Introduction 

Cross sections for the ionization of highly charged heavy 
ions by light neutral atoms are necessary for the calculation 
of relative abundances of various charge states in low energy cosmic 
rays. Theoretical values for these cross sections and an assessment 
of their reliability in the 0.1-100 MeV/nucleon regime are particularly 
important because of the paucity of experimental values for some 
of the cases of astrophysical interest. For these reasons we have 
calculated single electron loss cross sections in the first Born 
approximation (FBA) for hydrogen-like and helium like ions, as well 
as for ions with more than two electrons. In our calculations, the 
ejected electron is described by screened hydrogenic wavefunctions 
arid the neutral target atoms are characterized by Hartree-Fock form 
factors. 

The present application of the FBA to ion-atom collisions follows 

1 2 

closely the methods developed extensively by Bates and his co-workers. ’ 

1 -7 

Most of the earlier theoretical work using the FBA to calculate 

ionization processes was applied to hydrogen-like and helium-like 

systems being ionized by hydrogen and helium atoms. There have 

fi R-l 1 

also been several FBA calculations ’ for electron loss by 
H and He in heavy neutral targets, which treat the target using 
a closure approximation for the infinite sum over final target 
states. However, few calculations are available for heavy projectiles 
with more than two electrons colliding with neutral atoms. As a 
result, scaled binary encounter approximation (BEA) ionization 


cross sections have often been used to calculate relative abundances 

of the ionic charge states in beams of low energy cosmic rays. 

Scaled BEA cross sections have the advantages of availability in 

the literature and ease of calculation?, however, to use this 

approximation, one must argue that the structure of the target 

atom, usually H or He, is not important in the ionization process. 

The validity of such an argument is not completely clear in all 

cases. We have therefore included a comparison of the present 

12 

FBA results with those from the scaled BEA to help determine the 
latter method's reliability and accuracy. 

We also have examined the scaling properties of the cross 

sections for hydrogen-like ions with respect to both target and 

. +25 

projectile nuclear charges for projectiles ranging from H to Fe , 

and for neutral targets ranging from H to Fe. Finally, we have com- 
pare: our results for the ionization of N +1 (i=0-6) and 0 +l (i=3-7) 
by He with experimental data in the energy range from 0.1 to 10 
MeV/nucleon. These comparisons have proved very useful in deter- 
mining the expected region of validity of the FBA and the BEA in 
cases necessary fur the study of low energy cosmic ray ions, for 
which no data is available. 


II. Method of Calculation 

In our calculations, we assume that the major contribution to 
the loss of a single electron is direct Coulomb ionization. We then 
write the total ionization cross section as a sum over partial cross 
sections for each occupied subshell of the projectile; 

wvi-e f «JM' <" 

where N is the principle quantum number of the highest occupied shell, and 

T 

V is the relative ’'elocity of the projecti le- target system. Each a n ^ can 
be separated into two parts, one in which the target remains in the ground 
state, L and one which is the sum over all inelastic target processes, 

°ni- as 



( 2 ) 


After Fourier transforming the potential in the usual expression for the 

FBA cross section, one obtains for (V) (see, for example refs. 2,4,10); 

n 


and 




(V)= ?ira; 



-k 


a,X 



In the above expressions we have: a Q the Bohr radius, v Q = ac the Bohr 

velocity, Zj the target nucleus, and q = , 
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where ty and are the initial and final momenta in the 
center of mass system in a Q 1 . 

The target elastic form factor F(q) is given by 




( 5 ) 


(summed over occupied spin-orbitals), where. 


f:./^ si- <x(r)\e r ~\rjin> r , w 

in which ^(r) is a single particle spin-orbital for a target 

electron with coordinate r measured from the target nucleus. The 
incoherent scattering form factor, S (q ) in (4), is defined by 



where N-y is the number of target electrons (Ny - Z-y for neutrals). 

The function 8 n ^ (k,q) for the ionized projectile is given by^ 

0 ) kKZple ^njm^ <«) 


which is just the inelastic form factor for a bound to continuum 

transition, averaged over m, integrated over the angle k of the 

ionized electron, and multiplied by the number of electrons in 

the subshell, N n The wave functions are hydrogenic functions 

i f 

with effective charges Z p for the initial state and Zp for the final 
state. The states } T< , Z > and |n£m,Zp > have been constructed to be 

r r 
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orthogonal^ 4 in order to avoid additional terms in (8) from the 

interaction of the projectile nucleus with both the target nucleus 

f i 

and the target's electrons, which would otherwise appear when Zp f Z^. 
The effective charges can be specified for each pair of quantum 
numbers (n,£). 

We note that if Z^ = ll, then *(k,q) can be obtained by 

p p 

4 i 

defining Qe q/zj and using ^(k ,Q) evaluated with zj = zj = 1 in (8) 
above. The eqs. (3) and (4) can be rewritten in terms of Q. From 
this it is seen 4 that in terms of the scaled velocity, • 


V s V/ (Zp v 0 ) , 


(9) 


we find the following approximate relations for v»l: 


and 



oC 




(v) OC 


, . L 2 ,2 

Z T /(Z P v) . 


( 10 ) 


( 11 ) 


Finally, we discuss the limits of integration appearing in (3) and 
(4). The momentum q takes on its maximum and minimum values, q ^ and q-j , 
when k = 0: 


^x =K i ±K i[ 1 '^Kf 
^K l ± K l [: 


, -A 


ri p +Ag] i/a 

(Ip + A l t )] , 


(12a) 

(12b) 


6 


I 
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with ^ being the reduced mass of the projectile-target system, 

m the electron's mass, I the ionization potential in Rydbergc 
6 P 

for the (n,£) subshell of the projectile, and aEj the change in 
the target's internal energy. For q-, and q 2 in (3) aE. t = 0 and 
we have 

C| i ^ Ai I p ~ Ip/(fi V Z p ) (13) 

and 

^ ~ 2Kj_ . (14) 

For q^ and q 2 in (4), we must use an effective target excitation 
energy for aE 7 » since the expression for o* n(V) was derived by 
using closure to sum over all inelastic target states. We have 

g 

adopted the procedure of Lodge and let 

A E ® I + A 2 . 05) 

I r* 

Alternative choices for aEj have been studied previously^; also, 

1 5 

corrections to the closure approximation have been calculated. 

Using (15) in (12b) we obtain 

~(I p + I T )/(sv £p) 06) 

In the actual calculation we take q 2 , qg + “• The upper limit of the 
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k integration, k_ v in (3), is given by 

MlaX 

r 0 -, 1/2 

^ = )- I pj > ' 17 > 

while (15) gives 

, _l/2 

V x =[^(^Ki-^)-KVi T )J 081 

For k; ax In (4). 


,.t 



8 



III. Results 

In tills section we will describe tki results of applying 

the theory outlined in the previous section. In order to carry out 

1 fi 

the calculations, we have used analytic fits to Hartree-Fock 
elastic form factors for F(q) in (3). For the incoherent scattering 
factor defined in (7), we have used those obtained by Cromer. 17 
The projectile inelastic form factor in (8) was obtained as in ref. 

(14), using screened hydrogenlc wave functions. For each subshell 
we have used effective charges Z^ calculated by Block and 
Mendelsohn. Other choices of effective charges ’ would alter 
our results by an amount which is smaller than the error associated 
with the experimental data in most cases considered here. The 
ionization potentials for each subshell were obtained from Moore's 
Tables. 21 

A. Hydrogen-like Projectiles 

Fig. 1 contains the results of calculations for electron loss 

from H*" 1 , 0 +7 , and Fe +2 ^ in collisions with both He and C targets. 

To examine scaling with respect to Z p (see eqs. (10) and (11).) we have 

plotted Z^=Zp times the total cross section in cm 2 versus the scaled 

velocity, v (eq. (9)). The H-He data shown for comparison is from 
22 23 

Toburen et. al . and Stier and Barnett . These He results are 

O 

similar to those obtained by Dmitriev et.al. ('iee fig. 3b of ref. 3). 
The carbon target results have been displaced to the right by one decade 
(upper scale) for clarity. The H-C data was deduced from data for 
H on Hg, Og* COg, CH^, C 2 H 2 , CgHg, and using the sum rule for 

the measured cross sections, neglecting molecular effects. 
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We see that Z scaling for projectiles of higher charge is quite 

r 

good for He targets, but for C targets the scaling is reduced. The 
scaling for high Z p can be understood from the fact that q^ and qj of 
eqs. (13) and (16) are proportional to Zjr. Thus the minimum momentum 

r 

transfer required for ionization is much larger for high Z p , and 
the q integration of eqs. (3) and (4) covers a range in which only 
the tails of F(q) «1 and S(q) = 1 are seen. Thus, for large Z p , 
scaling is quite good for He, while for C, the influence of F(q) 
and S(q) extends to large q values, reducing the scaling. 

Comparison of the H-He and H-C calculations with experiment 
indicates that the FBA for heavier targets is inadequate near the 
peak of the cross section (v = 1). Similar results have' been 
obtained by others 8, 9,11 . Walters 10 has treated this problem in 
terms of an exact calculation of ct-j for ionization in the static 
field of the target. The inelastic target contribution o^ 1 was 
still calculated in the FBA. His results for H-A agree much 

E 

better with experiment than do results based on the FBA for cr-| S , 
however this method requires considerably more computational effort. 
In Fig. 2 we show the results of ionization cross sections for 

+25 

Fe losing its electron as a result of colliding with H, He, C, 
and Fe target atoms. The dashed line is the binary encounter 
approximation (BEA) as tabulated by Hansen, which we have 
extrapolated beyond the range he has given. The dashed-dot line 
is the FBA result for protons as targets (FBAP), which scales 

9 

exactly as Zy for other bare nuclei as targets. The ordinate is 
4 2 

(Zp/Zy) times the total ionization cross section, and the abscissa 
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is v, eq, (9). Thus this set of curves tests the scaling of the 

total cross section, cJ eq. (2), with respect to the target charge. 

2 El 

Actually the Zy scaling is appropriate for a , while a scales as 

Zy(c.f. eqs. (10), (11)), The BEA and F3AP curves scale exactly 
2 

as Zy, since bare nuclei don't yield a contribution to the total 
cross section corresponding to cr 1 for the neutral targets. Although 
the latter curves are not truly comparable to results for ionization 
by neutral atoms, we have included them for the reason discussed in 
Section I. As seen in Fig. 2, for the Fe +25 ionizing in H, the BEA 
is ^ 40% lower than our results for v = 1. At higher energies the 
discrepancy increases. 

-2 E II 

We have found that Zy a and Zy a obey scaling quite well 
if scaled separably. We also find that the curve for protons as 
targets? FBAP, is nearly identical to the elastic target cross section, 
o ,for neutral hydrogen below v s 3; for v = 20, the FBAP result 
exceeds a by ^ 47%. This result is expected, since, for charged 
targets, the long range Coulomb force gives rise to an E" 1 log E 
energy dependence, in contrast to the E dependence for neutral 
targets. The fact that the lighter targets give scaled ionization 
cross sections which are above the FBAP curve, while the heavier 

targets give values below the FBAP curve is explained by the fact that 

IE T T 

a / a “1/Zy. Thus the relative contribution to 0 from 0 , eq. (2), 

is reduced for heavy neutral targets. 


11 


B, Multi -electron projectiles 

In Fig. 3 we give the FBA results for total ionization 

cross sections as a function of energy calculated for several 

charge states of oxygen being ionized by He. Also shown are 

+2 

SEA cross sections for He targets and the data of MacDonald 

pc pc 

and Martin and Dmitriev et. al. ’ No error bars were 
displayed in ref. 23. However we have estimated the bars from 
the discussion given there and included them in our figure 
in order to aid in comparing our results to the BEA and to 
experiment, as well as to compare the four lowest energy data 
points of ref. (24) to the higher energy data. 

The FBA results (solid curves) were calculated by summing 

T f * 

the a p 's with respect to n and £, as in eq. (1), using Z = Z' 

n yC P P 

(see eq. (8)) in each case. For the lower charge states in which 

T 

the 2s electrons were involved, we did not recalculate a-j s with al- 

T 

tered screening and ionization potential since o-| s was only a 

small contribution to the total cross section as compared to . 

The effect of the additional screening by an outer electron can be 

+3 

seen in the case of 0 , where there is one electron in a 2p state. 

For this case we show both the result of the uncorrected sum, 

°ls + a 2s + °2p’ ( s0 ^ d curvs ) and ^ the sum in which o 2 g was 
recalculated (dash-dot curve) with the additional screening and 
altered ionization potential coming from the outer 2p electron. 

We conclude from this that such additional screening is only 
important for the lower charge states of the ion, in which the 2p 
subshell begins to fill. For heavier, more highly charged projectiles 
such as Fe, the effect of additional screening of the inner electrons 
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by the outer ones should be even less than in the case of oxygen. 

In comparing our results with the BEA in Fig. 3, we see that, 

+7 

with the exception of the 0 case, where the experimental error 
is quite large, our FBA calculation seems to be in somewhat closer 

•j*A +3 

agreement with experiment, especially for the cases of 0 and 0 . 

Also, there does not seem to be a systematic relationship between 

the BEA and our results for the difference charge states. 

In Fig. 4 we give the results of the FBA calculation along 

24 2R 

with experimental results of Dmitriev et. al. ’ for all the charge 

states of nitrogen in helium. The ionization cross sections were 

calculated in the same way that was described for the O'^-He collisions, 

2 12 3 

with the exception that for the electron loss from the 2s and 2p ’ ’ 

•f 

levels, l f l' (see eq. (8)). For these cases it was found that closer 

r r 

agreement with experiment was obtained if was chosen, as usual, to 
be the screeened charge for a given level, but Z was taken to 
be the asymptotic charge seen by the Ionized electron. Thus to cal- 

tL i J. 

culate the ionization of N , we took Z ' ? i (i + 1), for i <. 4. For 

P 

ionization of ions containing 2p electrons, we also recalculated the 
contribution from the 2s subshell with the additional screening from 
che outer electrons and the altered ionization potential. As in the 
Q-He cases, the Is subshell's contribution was not recalculated since 
I- was much smaller than the 2s and 2p contributions. 

From Fig. 4 we see that near the peak, the agreement of the cal- 
culated and experimental cross sections is quite good, especially for 
the N°, N + ^ , and N + ^ cases, ^or the N + ^ and N + ^ cases, the agreement 
at energies below the peak is rather poor; however, for energies near 
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the peak and higher, the agreement is quite satisfactory. In 

+6 +6 

view of the accord of theory and experiment for N and 0 , we 

found the discrepancy between our results and the data of ref. 24 
+5 

in the case of N somewhat surprising; however, this seems to be 

resolved by the results of a subsequent, more refined experiment and 

analysis. The triangles in Fig, 4 represent the "most probable" 

values of t 'n total ionization cross section determined by Dmitriev 

et.al. after considering the effect of metastable {Is, 2s) * S 

states remaining in the beam which reaches the collision chamber. 

The experimentally measured cross section, they find, is in general 

roughly a factor of two greater than the true cross section for 

ionization of helium like ions in their ground state. In the 

+5 

particular case of N this ratio appears to bring experiment and 
our F8A results into accord for energies near the peak of the cross 
section. 

+i 

As mentioned above, in order to fit the data for N , we had to 
take Z = + (i+1) for i _< 4. Fig. 5 illustrates the effect of 

r 

using the asymptotic charge for Z £ in the continuum state of ^(k,q), 
eq. (8), by comparing with the case in which Z^ = Z^, where 

j I Q 

Zp is the effective charge appropriate to the initial bound state. 

We have plotted total cross sections, (Z^, zM, for ionization 

n j*'- p p 

from the 2p^ and 2s^ subshells of N + ^ and N + ^, respectively. The 

T 

cross section o^p (2.0, 3.80) approaches a value ^ 40% greater than 
a 2p (3.80,3.80) at high energies, while a 2s (4.0,5.05) approaches 
a value ^ 20% higher than 0^(5.05,5.05) at high energies. Thus 
we see that, for ( Z^ , Z^), the difference between the two 

methods for choosing Z is significantly larger than the experimental 

r 
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error of * 20%, with (2.0,3.80) giving results in excellent 

+1 

agreement with the data for N in Fig. 4. 
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IV. Discussion 

In order to summarize our results and to put them into 

perspective, we have adapted a figure from the review article 

by Madison and Merzbacher. Fig. 6 is a schematic representation 

« * 

of the regions in velocity (v =V/(zJ v Q )} and charge (Zy/zJ ) "space" 
in which the plane wave Born approximation (PWBAeFBA), semi-classical 
approximation (SCA), and molecular orbital approach (MO) are expected 
to be appropriate. Although this representation was originally designed 
for completely stripped ions colliding inelastically with neutral 
atoms, and although the various boundaries are, of course, not as 
well defined as in the figure, we have plotted those regions in 
this space which correspond to our calculations. The low velocity 
end of each of the lines and areas labeled by projectile-target 
pairs corresponds to the lowest velocity at which there is qualitative 
agreement between our results and experiment. Thus for H-C and 
H-He systems, the FBA cross section for ionizing hydrogen agrees with the 
data at velocities which are consistent with the usual criteria for 
the validity of the FBA: 

Z T /Zj «1 and Zy/zJ «v 

For the 0 +1 and N +1 projectiles on helium, we find that the accord With 
experiment extends to velocities lower than expected and the extent to 
which the PWBA region of validity overlaps that of the MO and SCA 
is apparent. This comparison of calculated and experimental results 
is important in establishing the regions of validity of the FBA for 
application to collision systems of astrophysical interest for which 


16 


little or no data is available. 

From a consideration of the results for the speci 1c multi - 
electron projectile-target systems that we have calculated, we con- 
clude that-for the ionization of highly ionized heavy particles by 
light atoms, the FBA should give very reliable values of cross 
sections for velocities corresponding to the peak of the cross 
section and higher. 
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FIGURE CAPTIONS 


Fig. 1. Projectile charge scaling of total ionization cross sections 

as a function of scaled relative velocity {eq. (8)) for hydrogen- 
like projectiles colliding with neutral He and C target atoms. 

( — ) (lower scale) the - scaled FBA result for electron 
+7 +?5 

loss by H, 0 , and Fe projectiles of nuclear charge Z = Z p 

in He. ( ) (upper scale) the same processes in C. (a) data 

from ref. 22, and (b) from ref. 23. 

Fig. 2. Total electron loss cross sections for Fe scaled as (Z /Zj ) 
for various targets as a function of scaled relative velocity. 

Here ZsZ = 26. ( ) present FBA results for neutral target 

atoms of nuclear charge Zp ( — ) the BEA values from ref. 13, 

and ( — ’ ) the FBAP result, both for bare nuclei as targets. 

Fig. 3. Single electron loss cross sections for the ions 0 +1 (1=3-7) in 

collision with neutral He as a function of laboratory energy (E/A) 

in MeV/nucleon. ( ) and ( ‘ ) present FBA results, ( — ) the 

BEA results for He +2 targets. Data from (a) ref. 25, (b) ref. 24, (c) 
ref. 26, containing corrections for metastable states in the 
beam (see discussion of Fig. 4 in the text). 

Fig. 4. Total ionization cross sections for the loss of one electron 
by N (i=0-6) in He, versus laboratory energy in MeV/nucleon, 

{ ) the present FBA results (a) data of ref. 25; (b) from ref. 

26, after correction for metastables in the beam. 

T f 1 

Fig. 5. Comparison of total ionization cross sections, c* n ^ (Z p , Z p ), 

■j* «£ ^ 

resulting from choosing either Z p s Z p or Z p - (i+1) for the 
(n,£) subshell of N +1 (i=l,3). a) 2s 2 subshell of N +3 ; ( ’ ‘ — ) 
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a 2s T (5.05, 5.05), ( ) a gj (4.0, 5.05). b) 2P 2 subshell 

of N +1 ; ( — ) a 2p T {3.80, 3.80), ( ) a 2p T (2.0, 3.80). 

Fig. 6. Schematic of the regions in ciiarge and velocity space for which 
the molecular orbital approach (MO), the plane wave Born 
approximation (PWBA), and the semi -classical approximation (SCA) 
are expected to be applicable. The regions labeled by the 
projectile-target pairs N +1 '-He, 0 +1 '-He, H-C and H-He are for 
the results of the present work which are in accord with 
experiment. 
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